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ABSTRACT 



We test competing models that aim at explaining the nature of stars in the Milky 
Way that are well away ( \z\ ■> 1 kpc) from the midplane, the so-called thick disk: the 
stars may have gotten there through orbital migration, t hrough satellite me rgers and 
accretion or through heating of pre-existing thin disk stars. Sales et al. ( 20091 ) proposed 
the eccentricity distribution of thick disk stars as a diagnostic to differentiate between 
these mechanisms. Drawing on SDSS DR7, we have assembled a sample of 34,223 G- 
dwarfs with 6-D phase-space information and metallicities, and have derived orbital 
eccentricities for them. Comparing the resulting eccentricity distributions, p(e\z), with 
the models, we find that: a) the observed p(e\z) is inconsistent with that predicted 
by orbital migration only, as there are more observed stars of high and of very low 
eccentricity; b) scenarios where the thick disk is made predominantly through abrupt 
heating of a pre-existing thin disk are also inconsistent, as they predict more high- 
eccentricity stars than observed; c) the observed p(e\z) fits well with a "gas-rich merger" 
scenario, where most thick disk stars were born from unsettled gas in situ. 



Subject headings: Galaxy: evolution - Galaxy: kinematics and dynamics 
structure - Galaxies: individual (Milky Way) 



Galaxy: 



1. Introduction 



The common presence of so-called 'thick' stellar disk components in spiral galaxies, including 
in our Milky Way, has been established for some decades (jBurstein I Il978l ; iGilmore Reid I Il983l ; 
Yoachim &i Dalcanton 1120 06): the vertical distribution of stellar brightness or mass is more ade- 
quately fit with a thin - thick disk superposition, rather than with a single (thin) disk component, 
where the thick disk has a scale height that is typically three times larger. Thick disk stars further 
differ from stellar populations closer to the galactic midplane in terms of age, metallicity and rota- 
tional velocity. Certainly in the Milky Way, thick disk stars are generally older, more metal-p oor, 



and rotate more slowly around the galactic center (e.g. iFeltzing &; Bensbv 



2008 



Reddv 112009 ). 
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There is no consensus on the origin of the thick disk. Qualitatively different scenarios for the 
origin of thick disks have been discussed and cast in simula tions, resulting in d iffering chemical and 



kinematic properties for the thick disk stellar population. ISales et al. I (|2009l ) examined examples 



of these different scenarios and proposed the orbital eccentricity distribution of stars at 1-3 thick 
disk scale heights from the plane as a model discriminant for the case of the Milky Way. Broadly 
speaking, two of the four scenarios in ISales et al. I (|2009l ) presume that the thick disk stars got 
heated from a once thinner disk, and two explain the thick disk as a consequence of material that 
was deposited in the course of a minor merger. 

The first of these scenarios, radial migration, is based entirely on internal processes, and the 
thick disk is cre ated from stars migrating outwards from the kin ematically hotter inner regions of 



the Milky Way (jRoskar et al 



2008 



Schonrich &: Binney 1 120091 ) . In the second scenario, heating, 



the thick disk is predominantly t he result of rapid heating o f a pre-existing thin disk through one 
relatively massive merger event (IVillalobos fc Helmi 2008). I n this scenario, a modest fraction 



of the original thin disk is preserved (jKazantzidis et al. 



2008 



Villalobos fc Helmi 112003 ). and a 



small portion of present-day thick disk stars were formed in the satellite galaxy and are on highly 
eccentric orbits. 

In the accretion scenario, thick disk stars mostly form in an external satellite galaxy that gets 
disrupted while on a prograde orbit near the disk plan e. This can produce many properties of 
observed old thick disk components (jAbadi et al. Il2003l ). Finally, the gas-rich merger scenario is 
a variant of this idea, where a minor merger epoch brings gas into the galaxy from which (thick 
disk) stars form before the gas completely settles into a thin disk. 



Analyzing four concrete models based on these scenarios, ISales et al. I (|2009l ) proposed orbital 
eccentricity as a comparative diagnostic. For each simulation, they made eccentricity distributions 
for simulated stars distant enough from the midplane to eliminate thin disk stars. They pointed 
out that the eccentricity distribution p(e\z) differed significantly between these four scenarios (re- 
produced as histograms in their Fig. 3). 

The aim of this Letter is to compare these predictions to the observed eccentricity distribution 
of thick disk stars in the Milky Way, drawing on G-dwarfs from SEGUE DR7 with 3D position 
and 3D velocity information. At face value, this comparison favors the gas-rich merger scenario as 
the most important mechanism of thick disk formation. We also point out that the relationship 
between height, metallicity and eccentricity of thick disk stars has great potential for more stringent 
data-model comparisons with existing data. 



2. Data 



Sales et al 



2.1. Sample Definition 
(|2009l ) calculated their diagnostic eccentricity distribution in "solar neighbourhood 
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volumes" , which they devised as cylindrical shells between two and three scale radii of the model 
galaxies' thin disk (2 < R/Rd < 3). To reject most of the thin disk stars, they further constrained 
the stars to be a t 1 < | Scaled I < 3 (-z sca i c d = z/zq). For the vertical thick disk scale height zq, 
Sales et al. I (|2009l ) used ~ lkpc, except for their accretion model, which had zq = 2.3 kpc. 



It is impossible to match such a volume with any set of actual observations. Here, we draw 
observations mostly from the conical volume (with the Sun at its apex) that is naturally provided 
by the SDSS North Galactic Cap surv ey area. Specifically, w e drew up a sample of 34,223 G- 



dwarf candidate stars from SDSS DR7 (jAbazaiian et al. 112009 ) . for which spectra were taken and 
spectral parameters had been devised. Our selection criteria are outlined in Table [TJ and we 
provide the reader with the full data in Table [2] (avail able electronically o nly) . G dwarfs are the 



best-sampled spectroscopic target category in SEGUE (jYanny et al. M2009I ) . and a good fraction of 



these spectroscopic targets falls into the distance range of 0.8-6 kpc that is of particular interest 
here. 

Photometric distances (and errors) for these targe t stars were derived based on their ugr colors 
and apparent magnitude s, following the prescr iption of llvezic et al. I (|2008l ) which we expect to yield 



systematic errors < 5% (jKlement et al. II2009I ). These heliocentric distances allow us to derive the 



distribution of heights above the plane for these stars. For the model comparison, we n eed a sample 
of sta rs whose distribution p(z) reasonably matches exp{— z/zq), with zq ~ 1 kpc (as in lSales et al. 



20091 ). We found that in the range 1 < |-z SC aied| < 3, the p(z) of our sample matches exp(— z/zo) to 
within 20%; apparently, the SDSS selection function, the G dwarf luminosity distribution and the 
conical volume compensate to provide us with the appropriate height distribution. 



2.2. Phase-space coordinates 



The directly observed quantities for each star are its position, photometry, proper motion 
and its spectroscopic parameters, with the line-of-sight velocity and the metallicity [Fe/H] most 
relevant here. These quantities need to be translated into the phase-space coordinates r and v 
and their errors in the Galactocen tric reference frame. This is done using conversion matrices 
from Ijohnson Sz Soderblom I ( 19871 ). with v c = 220 km.s -1 for the circular velocity of the Local 



Standard of Rest (LSR) and 8 kp c for the Sun's distanc e to th e Galactic Center. The errors 
for each Vj were derived following Uohnson &: Soderblom I (119871 ). using the radial velocity and 
proper motions provided in DR7 and distance measurements taken from above. Typical errors 
are ~ 5 km.s -1 for the line-of-sight component ( Yannv et al. 2009)0, and ~ 25 km.s -1 for each 
transverse component at a distance of 2 kpc. The errors for each galactocentric ?j are calculated 
simply from the distance error under the assumption that stellar galactic longitude and latitude are 



Our formal radial velocity errors obtained from the database are smaller by a factor of ~ 2; however, this only 
has a minor effect on the spatial velocity components compared to the proper motion errors. 
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Table 1. Selection criteria for SEGUE data. 



Parameter 


Criterion 


(g - r) 


0.48 < {g-r) < 0.55 


{u - g)o 


0.6 < (u-g)o < 2.0 


{r - i)o 


-0.1 < (r-i)o < 0.4 


Absolute velocities 


\v\ < 600 km/s for v x , v y , v z 


T eff ,[Fe/H],[a/Fe] 


have valid values 


log(g) 


log(g) > 3.75 


E(B - V) 


E(B - V) < 0.3 



Note. — The three first rows impose color restric- 
tions valid for G stars. The fourth removes stars that 
are gravitationally unbound. Row five specifies the 
need for T e ff,[Fe/H],[a/Fe] data. Row 6 imposes 
high g to select dwarfs for which photometric distances 
can be estimated. The bottom row minimizes interstel- 
lar extinction. 



Table 2. Observed and derived properties 
of our 34,223 SEGUE stars. 



Note. - Available online clS Si machine- 
readable table. 
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known precisely. The importance of both the individual uncertainties and of possible systematic 
uncertainties in the distance scale or the circular velocity at the solar radius are assessed through 
their impact on the derived eccentricity distribution in the next section. 



3. Analysis 

3.1. Eccentricity Estimates 

In order to calculate the orbital eccentricities for each star, we need to adopt a gravitational 
potential. Here we chose a simple logarithmic potential, $ = v% * lnr, where r is the distance from 
the Galactic Center to the star in spherical coordinates. With this potential, we can calculate for 
each star i its total angular momentum Lj, its energy Ei, and the effective potential: 

*eff(r i ) = *(r < ) + ^|. (1) 
2rf 

At the peri- and apocenter of the orbit of star i, the energy equals the effective potential: 

E% = $eff(»*apo/peri) > ( 2 ) 

where we solved for the two roots r ap0/ / peri in equation by simple bisection. The eccentricity is 
then defined and given as 

e = rap ° ~ rpcri . (3) 

?"apo T ?*peri 



The resulting distributions of observed orbital eccentricities are shown - as red histograms - 
in Fig. [H taken over height ranges 1 < | Scaled I < 3. 



3.2. Eccentricity Uncertainties 



Before discussing these distributions, we assess the impact of two important error sources on 
these distributions: the individual phase-space uncertainties and the choice of the potential. For 
each stellar (r, v)j with its associated (<5r, <5v)j, we created 100 Monte-Carlo realizations of phase 
space coordinates, calculated the eccentricities and created 100 realizations of the eccentricity 
histogram. These showed that the individual measurement uncertainties result in little variance in 
the eccentricity histogram and also cause no significant bias (e.g. away from small eccentricities). 
We also explored what the impact of a 5% systematic uncertainty in the distance (from the Sun) 
would be and found the effect on p(e\z) to be small. We also compared p{e\z) resulting from our 
choice of the simplistic logarithmic potential to the distribution produced by a more complex (and 
realistic) potential, consisting of a logarithmic halo, a Hernquist bulge and a Miyamoto-Nagai disk 
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(e.g. iJohnston et al. Ill999l ). The one-to-one correlation of the resulting eccentricities showed very 
little scatter for eccentricities less then 0.6. For higher eccentricities, the two estimates still show a 
good correlation, but with some net bias to larger eccentricities in the three-component potential. 
Adopting this more complex form for the potential would leave p(e\z) largely unaffected for e < 0.6 
and slightly boost the high eccentricity tail. We also explored the impact of changing v c by =t 10 
km.s -1 in the logarithmic potential context and found it to have no significant influence on the 
eccentricity distribution. This is because v c enters both into the transformation to the Galactic 
rest-frame and into the gravitational potential. 



3.3. Cutting the Sample 



To focus on the eccentricity distribution of thick disk stars and to minimize the contribution 
from the halo, we only consider stars with tangential velocity v$ greater than 50 km.s -1 in the 
Galactocentric rest-frame system (i.e. we remove all stars on slowly-rotating and retrograde orbits, 
consistent with lSales et al. M2009I ) . The effect of this cutoff on the overall shape of the eccentricity 
distribution is negligible. For the histograms in Fig. [TJ we also eliminated the small fraction of 
stars with metallicity [Fe/H] < —1.2, as they may be chemically attributable to the halo. This does 
not affect the overall result except for somewhat reducing the incidence of very high eccentricities 
(e > 0.6). We return to the question of the corr elation between kin ematics and metallicity at the 
end, though it was not explicitly considered by ISales et al. I (|2009i ). Together, these restrictions 
remove 10.6% of stars in the original data set. 



4. Results 

We now turn our attention to the main result of this Letter, the comparison of the observed 



Sales et al. 


( 


2009 


); this 


order as in 


Sales et al. 



(2009, their Fig. 3). The red histo grams show the o bserved p(e\z) for 



the four different scenarios, each scaled to zq (Table 1 from ISales et al. II2009I ). For all but the 
accretion scenario (top left panel) these distributions look very similar, with a peak at e ~ 0.25 
and a pronounced asymmetric tail towards high eccentricities, e ~ 0.9. Unsurprisingly, for the 
accretion scenario the mode of the distribution is shifted towards higher eccentricities, as we sample 
the kinematics from 2 kpc < z < 7 kpc from the disk midplane. 

Figure [1] shows that neither one of the two processes where the thick disk consists predominantly 
of stars that were "kicked up" from a once thinner configuration matches the observed p{e\z). 
Orbital migration (bottom left) predicts an eccentricity distribution that is too narrow. The data 
show a higher fraction of stars both on nearly circular and on highly eccentric orbits. The heating 



-7- 



Eccentricity distributions for 1 < z 3C < 3 
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Fig. 1. — Eccentricity distributions at a given height range above the midplane. 



scenario (top right), where much of the thick disk was puffed up from a thinner disk through a 
massive merger, provides a reasonable match to the observations for e < 0.6. However, this scenario 
also predicts that the stars from the satellite involved in this merger should be part of the thick 
disk, but at very high eccentricities 0.7 < e < 1; the data show now evidence for this. In the 
accretion scenario (top left) much of the thick disk consists of deposited satellite debris; as this 
panels shows, this predicts too high eccentricities even for the bulk of thick disk stars. 

The predictions of the gas-rich merger scenario (bottom left) match the observations well, 
at least qualitatively if not in the formal sense: there is a maximum in p(e\z) at 0.25 and a tail 
extending to e = 0.9. Recall that in this scenario both the pre-existing disk and the merging 
satellite are very gas-rich. Therefore, e.g. the "culprit" stars from the satellite which imply many 
high-eccentricity stars in the heating scenario (top right panel of Fig. [1]) are absent. 
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Discussion and Conclusions 



We have compared the predictions made by four thick disk formation simulations to eccentricity 
distributions at \z\ > 1 kpc resulting from SEGUE data. Eccentricity pr oves a valuable para meter as 
it allows to clearly differentiate and assess the four scenarios outlined bv lSales et al. I (|2009l ). Direct 
accretion of stars from satellites generally predicts orbits that are too eccentric. Radial migration 
fails to generate the observed number of very low and very high eccentricity stars. Kinematic 
heating of a primordial thin disk, on the other hand, predicts too many stars with highly eccentric 
orbits accreted from the disrupted satellite. The accretion of gas onto the disk as a result of satellite 
mergers is generally consistent with the observed distribution. 

Taken at face value, the comparison of the eccentricity distributions argues strongly against 
three of the four scenarios; they favor an origin of the thick disk {\z\ > 1 kpc) from a series of 
gas-rich (minor) mergers, where most thick disk stars formed in situ and a high eccentricity tail 
arises from accreted stars. However, the simulations underlying three of the four scenarios were 
aimed at elucidating the physical processes capable of producing a thick disk component in late- 



type spirals; they w ere not necessarily tuned to match Milky Way properties (jBrook et al. 112004 



Abadi et al. 1120031 ). Therefore, they may not represent the full range of possible p(e\z) within 
these scenarios. Also, these scenarios are of course not mutually exclusive. For example, orbital 
migration must be present at some level, irrespective of the minor-merger history of the Milky Way. 
The bottom left panel of Fig.[T]only argues against orbital migration as the sole or dominant process. 



Speculation about variants of the seemingly rejected scenarios that could match the data is 
beyond the scope of this Letter. However, we want to point out the wealth of chemo-kinematic 
information that is now available to further test the origin of the thick disk. In Fig. [2] we use the 
G-dwarf sample to illustrate how well the distribution of stellar orbits within the thick disk can be 
mapped as a function of height as well as metallicity: as expected, orbits become more eccentric as 
one moves away from the plane. At a given height above the plane, the most metal-poor stars are 
on the most eccentr ic orbits, while sta rs on nearly circular orbits have metallicities closer to that 
of the Sun (see also llvezic et al. II2008I ). 



Such information provides further tests of the scenarios. E.g., IVillalobos Helmi I (|2008l ) 



find that in their thin disk heating simulation, the fract ional numbe r of ac creted stellar particles 
increases with height above the galactic midplane. And lSales et al. I (|2009l ) argue that regardless 



of the specific formation mechanism considered, accreted stars are always associated with the high 
eccentricity end of the distribution. 

Ultimately, the observed distributions, as e.g. in Fig. [21 should be tested against further 
predictions made by simulations in order to substantiate or refute gas-rich satellite mergers as a 
plausible thick disk formation mechanism in the Milky Way. 
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Eccentricity distributions 




Fig. 2. — Eccentricity distributions for a range of heights and metallicities. 
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